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Abstract 



Context. Stellar companions of accreting neutron stars in ultra compact X-ray binaries (UCXBs) are hydrogen-deficient. Their helium or C/O 

accretion discs are strongly X-ray irradiated. Both the chemical composition and irradiation determine the disc stability with respect to thermal 

and viscous perturbations. At shorter periods, UCXBs are persistent, whereas longer-period systems are mostly transient. 

Aims. To understand this behaviour one has to derive the stability criteria for X-ray irradiated hydrogen-poor accretion discs. 

Methods. We use a modified and updated version of the Dubus et al. code describing time-dependent irradiated accretion discs around compact 

objects. 

Results. We obtained the relevant stability criteria and compared the results to observed properties of UCXBs 

Conclusions. Although the general trend in the stability behaviour of UCXBs is consistent with the prediction of the disc instability model, in 
a few cases the inconsistency of theoretical predictions with the system observed properties is weak enough to be attributed to observational 
and/or theoretical uncertainties. Two systems might require the presence of some amount of hydrogen in the donor star. 

Key words, accretion, accretion discs - instabilities - stars: binaries: close 



1. Introduction 

It is well known that accretion discs around compact objects 
are subject to a thermal-viscous instability at temperatures 
corresponding to partial ionization of hydrogen. The model 
based on this instability explains the main properties of 
dwar f nova and lo w-mass X-ray binary (LMXB) outbursts 
(see Lasotal 2001 , for a review). In the case of transient 
LMXBs (often called soft X-ray transients - SXTs) X-ray 
irradiation of the accretion disc plays a fundamental role 



in th e stability criteria an d outburst physics (Ivan Paradij s 



19961 iKing & Ritten 1 19981: iDubus et alJ ll999t 1200 lb . An 



analogous instability should be present in ultra-compact 
binaries (P OT b ~ 60 min) whose accretion discs are kn own 



(2000) and IMenou et alj 020021) . However, since they did not 
include effects of irradiation, they were not applicable to 
UCXBs in which irradiation is important ( Del ove & Bild sten 
20031: iNelemans & Jonkerl2006l but they are of course relevant 
to AM CVn stars). 

The formation and evolutionar y channels of UC XBs are 
still subject to discussion (see e.g. INelem ans 2007, and ref- 
erence therein). In the so-called white-dwarf channel a dou- 
ble white-dwarf or neutron-star - white-dwarf binary evolves 
due to angular momentum losses through emission of grav- 
itational radiation. Such systems evolve from very short to 
longer orbital periods with ever decreasing mass-transfer rat e 
JPeloye & Bildste nl2003r.lDeloye et al.l2007r.lNelemansl2007l) . 



to be hydrogen deficient (see e.g. Nelemans 2005, 2007, for 
recent reviews). Since some AM CVn binaries (in which 
matter is accreted onto a white dwarf) show dwarf-nova type 
outbursts and several Ultra Compact X-ray Binaries (UCXBs 
- bin aries with a neutron-star accreting primary) are transient 
(see INelemans et alj|2006t INelemans & Jonkeii l2006) it seems 
that the thermal-viscous instability is operating also in these 
systems. Stabilit y crite r ia for hydrogen-de ficient discs wer e 
first obtained b y Sm akl d 19831) and later b y Cannizzo (11984 
Tsugawa & Osakil d 1997b . lEl-Khourv & Wickramasin ghe 
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In the second channel a binary composed of non-degenerate 
helium star and a white-dwarf or a neutron-star shrinks first 
to a minimum period then evolves back to long periods 
with a strongly decreasing mass-transfer rate. Finally, the 
third scenario involves cataclysmic binaries or LMXBs with 
evolved companions. Such systems evolve towards shorter pe- 
riods and when mass-loss from the companion uncovers the 
helium core their evolu tion is similar to the previous case 
( Podsiadl owski et al J2002b . n all three cases the mass-transfer 
rate decreases with orbital period. On the other hand the crit- 
ical accretion rate determining the stability of a hot accre- 
tion disc increases with the orbital period. This is because in 
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a stationary disc, the temperature decreasing with radius and 
the critical temperature being practically constant, the stability 
is determined by the physical parameters of the outer disc's 
ring. The disc radius in turn increases with orbital period. 
Therefore one should expect the shorter period UCXBs to be 
stable and persistent, while longer period systems should be 
transient. This expectation is roughly confirmed by observa- 
tions. However, the absence of stability criterion for irradiated 
hydrogen-deficient accretion discs has not allowed a detailed 
comparison of the evolution models and t he Disc Instability 



Mode l (DIM) with observations of UCXBs (Del ove & Bildsten 
20031) . 



In this article we study the stability properties of irradiated 
hydrogen-poor discs and apply our results to UCXBs. In sec- 
tion [2] we present stability criteria for helium discs and in the 
case of no irradiation compare our results with those obtained 
by other authors. Application to observed UCXBs is presented 
in section [3] and problems are discussed in the last section [4] 
In the Appendix we included updated results for hydrogen-rich 
discs. 



2. Thermal-viscous stability of helium discs 

We used an u pdated and slig h tly modified version of the code 
described in lHameurv et all dl998l) . iDubus et all {l999) and 



Dubus etal. (2001). 



2.1. Opacities and EOS 

The updated opaciti es were taken from OPAL 



(lgle sias&Rogerslll996l) and Ferguson et all d2005l) . OPAL 



Rosseland mean opacities run from log T = 3.75 to 8 and the 
Ferguson et al. opacities run from log T = 2.7 to 4.5. An aver- 
age of log k is taken in the temperature regime where the tables 
overlap, weighted by log T, so that the two tables connect 
smoothly. In addition to the opacities, the thermodynamical 
quantities for the gas m ixture (pressure, internal energy, etc.) 
are calculated following lPaczvriski dl969l) . Our accretion discs 
are pure helium (Y — 1) except where noted. In the cases 
wher e the metal composition was non-zero (Z + 0), we assume 
the Grevesse & Noels ( 19931) metal relative distribution. 



2.2. S-curves 

The stability condition is found by investigating the local ver- 
tical struc ture of the disc. The ve rtica l structure was calcu lated 
following lHameurv et alJ d 19981) and lDubus etaD dl999l) . The 



input parameters are Rw, T c , Ti n - and a (where R = Rw 10 10 cm, 
T c the midplane temperature, T m the irradiation temperature 
and a the viscosity parameter). The returned variables are the 
surface density X and the effective temperature T e s for which 
local thermal equilibrium is achieved. The solutions form the 
well-known S-curves in the (X, T e $, or T c , or M) plane. We 
computed several thousands of such S-curves for a between 
10~ 4 and 1, R between 10 6 and 10 11 cm, T m between and 
25,000 K, T c between 10 3 and 10 6 K. Example S-curves are 
shown in Fig.Q] The critical values of S and T e s are given by 
the two inflexion points in the S-curve. Their dependence on 
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Figure 1. Example S-curves for a pure Helium disk with vary- 
ing irradiation temperature 7V. The plot shows the effective 
temperature as a function of the surface density for a = 0.16. 
The different sets of S-curves correspond to different radii 
R = 10 6 , 10 9 and 10 10 cm. For each radius, the irradiation tem- 
perature TV is OK, 10 000 K and 20 000 K. The instability 
disappears for high irradiation temperatures. 

the input parameters can then be approximated by a numeri- 
cal fit, g iving the critical £. T c , J e ff as a function of a, R and 



as in 



Dubus et al. (2001). The critical mass accretion rate is 



derived from the above using 
8nR 3 



M = o-T 



eff 



3GM 



(1) 



The dependence on M is derived from the dependence on R 
since the S-curve is unchanged for a constant M/R 3 ratio. 

2.3. Non-irradiated helium disc: comparison with other 
calculations and impact of metallicity 

For the non-irradiated case T\ n — K the critical values for a 
pure helium disc (Y — 1) are 

T c + = 76000 a^R^Mf-^K 
T+ = 13100 ^ 01 ^ 08 M? 03 K 



ff _ 1J,UUU 0.1 "10 J "l 

M c + rit = 1.05 x 10 17 a ( ; 05 R 2 W 69 MT 0M) g s" 1 
2- rit = 1770 ff0 r^: 2 °Mr - 40 gcm- 2 
T- = 16000 tfo ! 14 R-£ 05 M^Mk 
r eff = 9700 ffo ( J 00 < 09 < 03 K 
M-, = 3.18 x lO^ao^^M^gs -1 



(2) 
(3) 



MO 



where Mi is the mass of the compact object in solar units and 
a = 0.1 ao.i- The superscript (+) designates a critical value for 
hot ("upper branch") discs, whereas (-) corresponds to the cold 
("lower branch") discs. 



other authors. 


Menou et al. 


Hameurv et al. 


1 19981 code) 


(7=1) 



|ive Ml-, for a pure-helium disc 



Ml, = 5.9 x 10 16 a'™ 1 R™ 2 M^ 87 g s" 1 . 



MO 



(4) 
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The fit to the critical mass accretion rate is very close to what 
we find except for the stronger dependence on a in the fit of 



Men ou et al.l (2002). This is the result of their assumpti on of 



a perfect gas equation of state: as discussed by Menou et al 



(2002), the dependence on a disappears when convective trans- 
port is included, as we have done. 

El-Khourv & Wickramas inghe 00) del the emission 



from AM CVn stars and present a few S-curves for some spe- 
cific values of a, R and M assuming a low hydrogen content. 
We find good agreement of the critical T e ff derived from these 
S-curves (see their Fig. 19). However, we find our critical S are 
hig her by « 50% . 

ISmaid(ll983l) who took Y = 0.98, Z = 0.02 obtains IT.. « 



790 g crrT 2 but in this case a non-stan dard way of calculating 
the cold branch was used dSmak 1984 , and private communi- 
cation). 



As noticed by Tsug awa & Os aki ( 1997), the critical values 
depend on metallicity (see their Fig. 3). This mainly changes 
the critical S, especially S~j . The reason is that at E~\ the 
opacity r ~ 1 and that changing the metallicity has a strong 
effect on the opacities at low temperature. It is therefore not 
surprising tha t our values for £~ . , fou nd using Y = 1, differ 
from those of Tsugawa & Osakif l 1997 ) whose baseline model 
uses Y = 0.97 Z = 0.03. Comparing the results of our numeri- 
cal fits (given above) to their S -curve for Z = 0.0004 (see their 
Fig. 3), we find that that our values of IT .. are within 15% of 

& n cnt 

their values. We also found very good agreement with their fits 
to the critical and S (Eqs. 4-7 of their paper) when using 
the same composition (Y = 0.97 Z = 0.03). 

One possible consequence of the metallicity dependence of 
IT ., is that running an unstable disc with a constant can lead 

cnt & 

to outbursts that have a larger amplitude when Z is low than 
when Z is high (because E~. /X+ is larger). This is intrigu- 
ing as large ratios of Elj./Eij. are known to be required t o ob- 
tain r ealistic lightcurves in the framework of the DIM dSmakl 
19841) . Typically, this is achieved by lowering a in quiescence 
i.e. assuming less efficient momentum transport in a cold disc. 
However, although the ratio of E~ VZ + ., increases with lower 

' ° cnt' cnt 

metallicities, the values are still far from the amplitudes re- 
quired for realistic lightcurves (about 10-20 compared to a ratio 
of about 3 here). Therefore this will not change the conclusion 
that one nee ds a lower a in quiescence for the DIM to work 
dSmaklll984h . 



2.4. Irradiated helium discs 

In det ermining the stabili ty of irradiated helium discs we fol- 
lowed |Dubiis_eialJ (11999b : we make a hypothesis on T[ n and 
find M for which the disc of given outer radius is stable. We 
assume that T^ n is given by 



o-Tl = C 



JAc 1 
AnR 2 



(5) 



where < C < 1 is the irradiation constant defined as in 
Dubu s et al. Jl999h . The value of C depends upon the radia- 
tive efficiency, irradiation geometry, albedo, irradiation spec- 
trum etc. 
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Figure 2. Example S-curves for a pure helium disc with vary- 
ing C. The plot shows the effective temperature as a function of 
the surface density for a = 0.16. The different sets of S-curves 
correspond to different radii R = 10 8 , 10 9 and 10 10 cm. For 
each radius, the irradiation constant C is changed from 10 , to 
10" 4 to 10~ 3 . Here, the irradiation temperature varies along the 
S-curve. Irradiation has little influence in the inner regions of 
the disc but greatly changes the S-curves at the outer radii. 

With Eq. (|5]l we calculate S -curves with as input parame- 
ters R\q, T c , C and a. Here, the irradiation temperature varies 
along the S-curve with the mass accretion rate required for ther- 
mal equilibrium. The critical point s are n ow given as functions 



of a, R and C as in iDubus et al. (1999). When one is inter- 



ested in stability (and not in time-dependent calculations) C 
is assumed to be constant for a given system. It is therefore 
more convenient to express the critical points as a function of 
C rather than of T m (which itself depends on R and M). We 
computed several thousands of S -curves for a between 10~ 4 
and 1, R between 10 6 and 10 11 cm, T c between 10 3 and 10 6 K, 
Mi between 0.1 and 100 M and C between 10" 10 and 10" 1 . 
Examples are given in Fig. [2] 

Irradiated pure helium disc 

For C < 10~ 6 one should use the non-irradiated fits given by 
Eq. |@J. For C > 10~ 6 the fits to the critical values are given to 
a good approximation (average relative error < 25% for Eij t , < 
10% for J e + ff and < 2% for logM c + it ) 

= 227C:° 16 < 79 < 96 -° 041ogC - 3 M- a25+0041ogC - 3 gcm- 2 



l 



r e + ff = 8730CI^ao i Ul R 



-0.06 „-0.19 D 0.02-0.011ogC- 3 ^0.02+0.01 log C-, K 
-uu.s ,-i)i)| ,_,-() 12-0.01 logC-j ^0.06+0.01 logC-j 



r c + = 4480c:rx.i *io 



10 



.0.20 „,-0.01 D 0.13-0.011ogC_ 3 ^-0.19+0.01 logC- 3 



K 
K 



T- m = 22600 C_ 3 a 0A R m 

HW+ --7 1 v in!6 r>"0.22 -0.03-0.01 log C_3 p 2.51-0.051ogC- 3 
m crit ~ X iU U -3 "0 1 K 10 

. M -0.74 + 0.051ogC_3 gs -l (6) 

where C = 10 3 C- 3 . 

We show only the results for the upper-branch (hot discs). 
Since in quiescent SXTs X-ray irradiation is negligible, the 
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relevant criterion for the lower-branch (cold discs) is the one 
without irradiation. For example, according to the DIM in a 
quies cent disc th e surface density must satisfy 2 < S~ (see 
e.g. lLasotalbOOlh . Since in such a disc M ~ R 265 (Eq. the 



ratio of viscous to irradiating fluxes F^/Fjry ~ R 165 and self- 



irradiation is never important (Du bus et al . 1999). 



Irradiated disc with mixed composition 

In order to discuss the influence of the hydrogen depletion we 
also calculated the case with X=0.1 and Y=0.9. We find 

V+ *A nr»-048 „-0.80 r,0.92-0.041oge- 3 ,,-0.2 1+0.05 log C- 3 

gem -2 

T; = 18700 Ct 07 «o°i 21 /?: o ao5 - ao21ogC - 3 M oo6 + o.o2iogc_3 K 

T+ /loen n-006 ^-0.02 D -0.16-0.021ogC- 3 ,,0.09+0.02 log C_ 3 v 
l eS = 4280 C_ 3 a Q1 K w M l K 

t+ nenn /-0.19 -0.02 D 0.09-0.021ogC- 3 ,»-0.16+0.021ogC-i v 
l in - 12500 C_ 3 a 01 rt 10 M l K 

M c + nt = 1.9 X 10 15 Clf 5 ^O^O.OllogC-3 ^2.38-0.061ogC-3 



M -0.64+0.071ogC_ 3gs -l 



(7) 



3. Application to ultra-compact binaries 



The neutron star's stellar companions in Ultra Compact X-ray 
Binaries (UCXBs) P 01D < 60 min cannot be hydrogen-rich stars 
since such stars would not fit into the orbit. The 10 confirmed 
UCXBs (see Table [3] have orbital periods in two ranges: one 
between ~ 10 and ~ 20 minutes, the second ~ 40 and ~ 50 
minutes. The four systems below 20 min. are classified as per- 
sistent X-ray sources (except for 4U 1543-624 they are X-ray 
bursters). Of the six binaries with periods above 40 min. one 
is classified as persistent: 4U 1916-05 which is an X-ray burst- 
ing dipper. The binary 4U 1626-67 (which contains a 7.7 s X- 
ray pulsar) used to be classified as persistent but it s brightness 
is slowly but clearly decaying dKrauss et al.ll2.007D . The other 
four systems are accretion-driven millisecond pulsars and are 
all transient systems. In general therefore UCXBs follow the 
expected pattern: shorter period are stable, at longer periods 
discs become unstable. However, as we shall see, comparing 
actual stability criteria with observations makes things a bit 
more complicated. In any case one should keep in mind that the 
division of UCXBs into persistent and transient might be "time 
dependent". We already mentioned 4U 162 6-67 (which should 
fade out into quiescence in 2-15 years, see iRrauss et al.ll2007l) 
but also the UCXB candidate 1H 1905+000 that was persistent 
for at least 1 1 years has now completely disappeared from the 
X-ray sky dJonker et al.ll2007l). 

According to IPaczvriski ( 1977 ) the maximum disc radius 
can be written as 
# D (max) 0.60 



a 1 +q' 

(valid for 0.03 < q < 1), where 



a = 2.28 X 10 y M} /3 (l + q) 

is the binary separation; P m ; n being the orbital period in min- 
utes. For convenience we will write the disc outer radius as 



cm 

min 



(8) 



(9) 



Rr 



2.28x 10Vm! /3 p 2 ( 3 cm 

J 1 min 




Figure 3. Stability crite ria for accretion discs in UCXBs. 
(Compare with Fig. 1 1 in Delove & Bildstenl2003l) . Four crite- 
ria are plotted (from top to bottom): (i) for non-irradiated pure- 
helium discs, (ii) for irradiated pure-helium discs, (iv) for an ir- 
radiated mixed-composition disc and (v) for an irradiated solar- 
composition disc. The accretion rates of persistent UCXBs are 
marked by circles. Transient systems are identified by triangles. 
Dashed lines correspond to three UCXB evolutionary tracks 
by Deloye (private communication). From bottom to top they 
represent evolution of binaries with a 1 .4M neutron star and 
a secondary with initial mass M2 and degeneracy parameter 
iff (see text) equal respectively: O.1M , 3.5; O.325M , 3.0 and 
O.325M , 1.5. 



where 



0.60 



(1 + 4) 2 / 3 ' 



(ID 



The value of C is not known a priori. Since it is a measure of the 
fraction of the X-ray luminosity that heats up the disc it con- 
tains information o n the irradia t ion ge ometry, X-ray albedo, X- 
ray spectrum, etc. iDubus et al.1 (119991) found that the observed 
optical magnitudes and stability properties of persistent low 
mass X-ray binaries were compatible with a value C ~ 5 x 10 -4 . 
It is not clear that the same value describes well the properties 
of UCXBs. Therefore will take C = 10~ 3 as the results do not 
vary much for the range 5 10~ 4 < C < 2 10~ 3 , say. With this 
value of C, from Eqs. ©, (JHJ and (0, assuming Mi = 1.4 
a = 0.1, / = 0.6, one obtains the following relation between 
the critical accretion rate and the orbital period for a pure he- 
lium X-ray irradiated disc 



M+ = 2.4 x 10- 12 /" 1 P^f'Moy" 1 . 

cnt -'0.6 min u J 



(12) 



In the non-irradiated pure-helium disc the equivalent relation 
reads 

M c + lt = 8.2 x 10- 12 / 2 fP m ™ M y- 1 . (13) 

For an irradiated mixed composition (X=0.1, Y=0.9) the criti- 
cal accretion rate is 



(10) Ml 



2.9xlO- 13 / 2 6 38 J P 1 ^M y- 1 . 



(14) 



Table 1. UCXBs (P orb < 60 min) 
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System 


P orb (rrun) 


M(M y-'^ 


Typel 


Comment^ 


411 1 820-30 1 


11.42 


5.1 x 10 9 


p 


(in GC) 


4U 1543-624 2 


18.2 


5.5 x 10~ 10 


p 




411 1850-087 3 


20.6 


1.3 x 10" 10 


p 


(in GC) 


M15 X-2 4 


22.58 


1.2 x 10" 10 


P 


(in GC) 


XTE J1807-294 5 


40.07 


1.9 x 10-' 2 


T 


(MSP) 


4U 1626-67 6 


41.4 


2.0 x 10- 10 


T 


(young pulsar) 


XTE J1751-305 7 


42.42 


4.5 x 10~ 12 


T 


(MSP) 


XTE J0929-314 8 


43.6 


4.1 x 10~ 12 


T 


(MSP) 


4U 1916-05 9 


49.48 


7.6 x 10" 10 


P 




SWIFT J1756.9-2508 10 


54.7 


9.3 x 10~ 13 


T 


(MSP) 



" see comments in the text 

6 P - persistent, T - transient. 

c GC - globular cluster, MSP - millisecond pulsar. 



i Zdziarski et al 



_ ll2007|) 2 dWang & Chakrabartvll2004l). 3 jsidoli et al.ll2006l). 4 jDieball et alj|2005l). 5 dMarkwardt et alj|2003l; iFalanga et alj 
l2QQ5h, 6 jKrauss et al j |2007|) , 7 dMarkwardt et alj 12001 iGierlinski & Poutanerj |2005|) , 8 jGallowav et alj |2002| ; ljuett et all l2003h , 9 
djuett & Ch akrabartv 2006), 10 dKrimm et al.ll2007h . 



and for the solar composition the critical accretion rates for 
irradiated discs is 



of he lium UCXB MSPs can be expressed as (see also lLasota 



<,t = 1-4x10- 



13 f 2.41 D 1.61 



■'( 



M y" 



These four criteria are plotted in Fig. [3] 



4. Discussion 



Deloye & Bildster] (2003) studied the structure and evolution 
of UCXBs and showed that for orbital periods > 30 min the 
mass-transfer rate cannot be > lO^M^y- 1 if they evolved 
adiabatically from systems filling their Roche-lobes at P or b ~ 
10 min. As seen in Fig. [3] one arrives at a similar conclu- 
sion using full stellar models (Deloye, private communication). 
These models are parameterized by the total binary mass and 
the initial degeneracy parameter of the secondary given by 
iff = Ep/kT c , wh ere Ep is the Fermi energy and T c the central 
temperature (see lDelove et al.ll2.Q07b . In general the evolution- 
ary tracks are consistent with the stability properties of indi- 
vidual systems but there are interesting cases where this can be 
questioned. 

All UCXBs known to be transient (the millisecond pulsars 
and the young pulsar binary 4U 1626-67) have mass transfer 
rates well below the stability limit for irradiated helium discs. 
However, a word of caution is needed as it is not clear that the 
criterion in form of Eq. ( fT2~b can be applied to MSP binary sys- 
tems as their sec ondary masses are supp osed to be < 0.0 1M , 
As mentioned in lYungelson et al.1 (12006b . for lower values of q 
matter transferred from the Roche lobe circularizes on unsta- 
ble orbits and it is unclear how disc formation proceeds. The 
actual outer disc radius may therefore be smaller (but not by 
much, see below) than we have assumed above, which would 
act to make the systems stabler than anticipated. On the other 
hand according to Eq. (fT2l the maximum outburst luminosity 



2008) 



(15) L n 



3.5 x 10 



37 



1.67 



erg s 



(16) 



which agrees well with observations. In fact in all UCXBs 
mass-ratios could be less than 0.03 but the fact that luminosities 
are consistent with the size of a standard disc model suggests 
that the stability criterion is applicable even in such extreme 
binary systems. 

4. 1 . Persistent systems 

The mass-transfer rate were calculated by using the luminosi- 
ties quoted in the references mentioned in Table [3] assuming 
a 1.4M neutron star with a 10km radius. Among the systems 
classified as persistent, two (4U 1820-30 and 4U 1543-624) are 
stable according to the stability criterion for irradiated helium 
discs but three other binaries (4U 1850-087, Ml 5 X-2 and 4U 
1916-05) should be transient according to this criterion but ap- 
parently are not. One can try to explain this apparent contradic- 
tion in three ways. 

First, the mass-transfer rates (bolometric luminosities) 
could be underestimated. Second, the outer disc radius has been 
overestimated. Third, the companion are not pure-helium stars 
but (still) contain some amount of hydrogen. Let us note that 
C/O - donors would make things even worse for the DIM as 
the corr esponding stability curves would be above th e helium 
curves dMenou et al. 2002 ; Delove & Bildsten 2003) because 
of their higher ionization potentials. 

Considering the first possibility one remarks that of the 
three outliers the two with P or b ~ 20 min are located in globu- 
lar clusters, so one could assume the distance to these sources 
is well known. However, the distance to the peculiar globular 
cluster NGC 6712 first believed t o be 6.8 kpc (ICudworthll988l) 
is now measured to be ~ 8 kpc bv lPaltrinieri et al.1 (1200 lb . Also 
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in the case of M15 X-2, the mass-transfer rate based on the X- 
ray luminosity used in Table 1 could be an underestimate and 



M > 4 x 1(T 10 M Q y" 1 dDieball et alj|2005l) . 



At first, the second possibility seems to be more promising 
since the dependence of stability criterion on the outer disc ra- 
dius is very strong. One could think that it is enough to take a 
value smaller than the maximum to get rid of the problem. For 
example at P 01 -b ~ 20 mn, an outer disc radius of 0.3 a (a factor 
~ 2 smaller than assumed above) is enough for the X-ray heat- 
ing of a helium disc to stabilize the system. However, already 
for q — 0.1 the circularization radius is equal to 0.3 a so for the 
likely smaller mass-ratio of systems of interest it is not a viable 
solution. 

It seems finally that the third possibility is the most likely 
way out of the stability problem for the two systems. A 
small fraction of hydrogen left (X > 0.05, see Fig. [3] l would 
make them stable. For 4U 1820-30 ICummind (120031) found 
that X ~ 0.1 is compatible with the X-ray burst properties 
of this system. This conclusion, however, does not apply di- 
rectly to the two systems of interest but evolutionary models of 
Podsiadlowski et al.l (120021) would allow them to have X > 0.1 
(rather paradoxically these models predict mass-transfer rates 
higher than observed, in fact in the helium stability range.) 

The evo l ution of 4U 1916-05 is rather controversial. 
Nelso n et al. suggested that mass-transfer rate deduced 

for 4U 1916-05 points to an evolved secondaries as donors 
in this system s. This would allow t he presence of hydro- 
gen. However, iNelemans et al.1 (|2006) challenge this conclu- 
sion finding a He donor and a high N abundance. However, it is 
not clear how this is consistent with the white-dwarf channel. 

Finally, the 20 min systems could have evolved thr ough 
the so-called "magnetic capture' ' (Ivan der Sluvs et al]2005i In 
such a (very unlikely) case they still contain enough hydrogen 
to fulfill the stability criterion for X < 0. 1 . 



The case of 4U 1626-67 is different. The neutron star in this 
UCXB is a young pulsar and it is transient in a different sense: 
its "outburst" does not last tens of days but tens of years and 
it is not clear at all that this behaviour has anything to do with 
the thermal-viscous instability of the DIM. The value of mass- 
transfer ra te given in Ta ble [3] is just a reference value given by 
Eq. (4) of iRrauss et al.l ( 12007b . The source has been on since 
1977 and seems to be decaying but the estimate that it will off 
in 2-15 years is based on dubious premises. One can say only 
two things: 

The mass-transfer rate of 2.0 x 10~ 10 M y" 1 i s inco mpati- 



ble with e volutionary track s of lDelove & Bildsten (2003). (See, 
however. lYungelson 2008 ). 

The mass accreted during ~ 30 years was ~ 5 x 10 24 g (for 
D = 3kpc) which is compatible with the maximum mass in a 
cold quiescent helium disc: 



I AM, 



(18) 



(see Eq. 3.) However, the recurrence time could be longer than 
1000 years, say and the mean accretion rate over the cycle 
< 2.5 x lO'^Mpy' 1 , c ompatibl e with evo l utiona ry tracks of 



Deloye & Bildsten (2003) and Deloye et al. (2007) 
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4.2. The unstable UCXBs 

Among the five systems in this category, four are accreting mil- 
lisecond pulsars (MSP). The mass-transfer rate was estimated 
from the formula 



M tran <4^(-) 2 J F 0Utb ^, 

V C I free 



(17) 



where D is the distance to the binary, F out b is the mean bolomet- 
ric flux during outburst; f ou tb and f rec are respectively the out- 
burst duration and the recurrence time. Except forXTE J1751- 
305 where the recurrence time is known to be 2-3 years, f rec was 
assumed to be 10 years and M tan i s therefore an uppe r limit. 
Our estimates our close to those of IWatts et al.l d2004l) . In all 
cases the distance to the sources is poorly known so in addition 
to the unknown recurrence times this makes mass-transfer esti- 
mates somewhat uncertain. From the point of view of the insta- 
bility criterion this is of no importance as all the MSP are com- 
fortably well below (two orders of magnitudes) the threshold. 
Fig.[3]shows also evolutionary tracks calculated by Deloye (pri- 
vate communication) and one mi ght suspect that MSP do nors 
are rather C/O than He stars (see lDeloye & Bildstenl l2003). 
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Appendix A: Solar-composition discs 

For comparison with results for helium discs we present here 
the fits to critical values for solar composition (X = 0.7, Y = 
0.28 and Z = 0.02) using the same procedure as described in 
the main text. 



Non-irradiated solar-composition disc 

= 39.9 a- o r°R\ U M^ gem- 2 

T: = 30000 18 sTAff- 01 K 

r e + ff = 6890/?7 009 M° 03 K 

M c + nt = 8.07 x 10 15 a ( ; 01 R 2 M M^ 9 g s" 1 (A. 1) 

^ lt = 74. 6ao r^!o 18 Mr - 40 gcm- 2 

t; = 824o4; 4 j R7 010 m; ,04 k 

r eff = 52io< 10 < 04K 

M c " rit = 2.64 x 10 15 a° f R 2 * s M^ 5 g s" 1 



Irradiated solar-composition discs 

The fits were obtained with the parameters varied in the range 
10~ 5 < a < 1, 10~ 4 < R w < 10, 0.1 < Mi < 100 and 10~ 6 < 
C < 1. The average relative uncertainty on log M^ rit is < 2%. 

v+ ot/o-0 28 -0.78+0.01 log C_3 D 0.92-0.071ogC_ 3 

. M -0.19 + 0.061ogC_3 gcm -2 

T: = 16300 CI°- 10 afi 11 R^- ^C- 3 ^.05+0.02 log c_ 3 K 

r e + ff = 404oc:r<^ro 015_00210sC "< 09+00210gC " K 

= 10500 C° 3 16 C^ 1 ^O-O-OllogC-, M -0.16 + 0.021ogC_3 K 

»>+ _ q c v i nl4 ^-0.36 0.04+0.01 log C-3 p 2.39-0.101ogC- 3 
M cnt ~ y J X iU U -3 "o.l 10 

_ M -0.64+0.081ogC- 3 gs -l 

(A.2) 



